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Nonuniform Spin Triplet Superconductivity due to Antisymmetric Spin-Orbit 
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We show that the nonuniform state (Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state) of the 
spin triplet superconductivity in noncentrosymmetric systems is stabihzed by antisymmetric 
spin-orbit couphng even if the magnetic field is absent. The transition temperature of the spin 
triplet superconductivity is reduced by the antisymmetric spin-orbit coupling in general. This 
pair breaking effect is shown to be similar to the Pauli pair breaking effect due to magnetic 
field for the spin singlet superconductivity, in which FFLO state is stabilized near the Pauli 
limit (or Chandrasekhar-Clogston limit) of external magnetic field. Since there are gapless 
excitations in nonuniform superconducting state, some physical quantities such as specific heat 
and penetration depth should obey the power low temperature-dependences. We discuss the 
possibility of the realization of nonuniform state in CePtaSi. 

KEYWORDS: CePtgSi, heavy fermion superconductivity, spin-orbit coupling, FFLO, spin triplet supercon- 
ductivity 
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1. Introduction 

Recently, superconductivity in the system without in- 
version center attracts much interest. The first heavy 
fermion superconductor without a center of inversion has 
been observed in CcPtaSi,^ and several pressure-induced 
superconductivity in the systems without inversion cen- 
ter have been reported, including CeRhSis,^ CelrSis^ 
and UIr.^° In CePtaSi antiferromagnetic order exists be- 
low Tjv = 2.2K followed by superconductivity below 
Tc = 0.75K. The upper critical field i?c2(~ 5T) exceeds 
the Pauli paramagnetic limiting field (or Chandrasekhar- 
Clogston limit), ^ which is thought to be an evidence of 
the spin triplet pairing. Recent experiments on Knight 
shift^' ^ also indicate the spin triplet superconductivity. 

On the other hand it has been believed that inversion 
symmetry is required for the spin triplet superconduc- 
tivity.'* Frigeri et al.^ have studied the spin triplet super- 
conductivity in the system without inversion symmetry 
in detail. They took the Rashba type spin-orbit coupling, 

(1) 



Hp ~ a 



gk • CTs,s'4,sCk,s' 
k,s,s' 



where cj^^ and Cks are creation and annihilation operators 
of electrons with momentum k and spin s, respectively, 
and gk = {—ky, k^, 0) as a model for CePtaSi, the point 
group of which is C4i, . They obtained that the spin triplet 
pairing is not affected by this spin-orbit coupling, if the 
d-vector of the spin triplet pairing satisfies d(k) || gk, 
while the transition temperature of the triplet supercon- 
ductivity is reduced if d(k) gk. 

Another interesting property of superconductivity in 
CePtaSi is the existence of the line nodes in the energy 
gap indicated by the temperature dependence of pene- 
tration depth^ and thermal conductivity.^ Since the spin 
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triplet state with d(k) || gk proposed by Frigeri et al.^ 
has points nodes at kx = ky = and no line nodes, this 
state is not consistent with the experiments. In order to 
explain the existence of the line nodes, the mixing state 
of triplet and singlet pairings has been studied. The 
strong spin-orbit coupling and the energy dependence of 
the density of states (particle-hole asymmetry) have to 
be assumed to obtain the mixed state of spin singlet and 
spin triplet. 

We propose another possibility, nonuniform spin 
triplet state, to explain the experiments in CePtaSi. We 
will show that there is a similarity between the pair- 
breaking effect due to the anisotropic spin-orbit coupling 
for the triplet superconductivity with d(k) H gk and the 
Pauli pair breaking effect due to the magnetic field for the 
spin singlet superconductivity. In the latter case Fulde 
and Ferrel,^^ and Larkin and Ovchinnikov^^ have shown 
that the Pauli pair breaking gives rise to a new pairing 
state between electrons (k-|- ^, t) and (— k-|- ^, |) on the 
exchange-split parts of the Fermi surface due to Zeeman 
effect of magnetic field, if the external magnetic field is 
close to the Pauli limit. This state is called FFLO state 
and it has been observed in CeCoIns,^'^ recently. In this 
paper we show that the similar pairing state to FFLO 
state is possible for the spin triplet superconductivity 
with d(k) 11 gk even if the magnetic field is absent. In 
the FFLO state without magnetic field the order param- 
eter is not uniform and the absolute value of the order 
parameter is zero in parallel planes with period in real 
space. Then some physical properties obey the power-low 
temperature dependences. 

2. Model and order parameters 

We adopt a single-band model with electron band en- 
ergy Ck measured relative to the Fermi energy, which is 
assumed to be spherical when there are no interactions. 
We use the weak coupling approach, taking the pairing 
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interaction as 



Full 



k,k',q s,s' 



(2) 



We assume that the interaction is finite and attractive 
close to the Fermi energy with the cutoff energy Cc and 
that the pairing interaction does not depend on the spin 
and the total momentum of the pair, q. This Hamiltonian 
satisfies time reversal and inversion symmetry. The ab- 
sence of inversion symmetry is introduced by spin-orbit 
coupling term, Hp (eq. (1)), with 



gk = 



3 1 



( ky,kx,0), 



(3) 



2kF 

as studied by Prigeri et al.^ For simplicity, we set 

(I gk |^)k = 1 where (• • • )k denotes an average over the 
Fermi surface. Then, the single-particle Hamiltonian is 



-f^o = ^ [^k + agk • d]ss' clicks'- 



(4) 



k,S,5' 



The normal state Green's function is obtained in the 
2x2 matrix form as 

G°{\s.,iuin) = G+(k,ia;„)(To -|- (gk • a)G_(k,za;„) 
Where 



(5) 



G±(k,ia;„) = 



1 



1 



1 



ek,± = ^k ± a I gk I, 
1 



gk 



gk 
|gk 



^k"^ + ky 



( ky,kx,0) 



(6) 
(7) 

(8) 



LUn = (2n + Vj-KksT is the Matsubara frequency, and cto 
is the 2x2 unit matrix. When the electrons (k + ^, si) 
and (— k+ 3, 52) make a pair, the linearized gap equation 
for FFLO state is written in a 2 x 2 matrix form^ as 

A,,,(k+ ^) = -UbT^Y, E ^1ck'G0,,(k' + f ,ia;„) 



k',n si,S2 

X A,,,,(k' + |)GO,,,(-k' + |,-ia;„). 

(9) 

The gap function is decomposed into a spin singlet part 
[V'(k + ^)] and a spin triplet part [d(k + ^)], as 

A(k + |) = [V'Ck + |)ao + d(k + |) . d]icTy. (10) 

Spin triplet part and spin singlet part are mixed in gen- 
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d(k+|) = -A;BTe^Fkk' 

n,k' 

([G'_G+g-G+G_g']*(k' + |) 
+ G+G+d(k' + |) 

-iG_G_(gxg')*(k' + |) 

+ i[G+G_i' + G_G+g]xd(k' + |) 

+ G_G_ [g x d(k' + |) x g' - g ■ d(k' + |)]) , 

(11) 

where we have used the short notation for the products: 



G„Gb = G„(k + f,io.„)Gfa(-k- 



-iujn) with a, & = ± 



and set g = gk'+f, g' = g_k'+a- We assume that the 
density of states is constant near the Fermi surface, i.e., 
the density of states is iV(0) even if the Fermi surface 
is shifted by the spin-orbit coupling. In this assumption, 
the particle-hole symmetry is satisfied and the singlet 
and the triplet order parameters do not mix. In other 
words, V'(k' -|- §) terms in the right hand side in eq. (11) 
vanish. Then we get the gap equation for triplet super- 
conductivity as 

d(k+^) 



= -ksTe E ^kk' (G+G+d(k' + |) 



n,k' 



+ i[G+G_g' + G_G+g]xd(k' + |) 



G_G_ 



g X d(k' + |) X g' - g ■ d(k' + |)g 



In the following we set q = (0, 0, kpq) and 

h^k%q cos, 6 



C±k-Fa = ^k ± 



2m 



(12) 



(13) 



8m 



-. Then we find the transition 



where ^k = & 
temperature for spin triplet pairing is given by 

m ric — l 



ln(-^) = -7rfcBTe(Im ^ 



n=0 



((/l_ + h- + /2+ + /l+ - 

- {.h- - /2- - /2+ + ,/l+) 



2(g'.d(k+|))(g.d(k- 



)) 



-g g 



d(k-Hf) 



where 



(14) 
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Basis function Order parameter 



A 


d(k) = 


/ 3 1 / ; ; n\ 




c{(\c\ - 
ai^Kj - 


- -^(k k k \ 




d(k) = 




A2 


d(k) = 




Bi 


d(k) = 




B2 


d(k) = 




E 


d(k) = 


= V3;^(fcz,0,0) 




d(k) = 





Table I. Order parameters in dy point group. The order param- 
eters are normalized as (|d(k)p) = 1 for the spherical Fermi 
surface. 

and Tco is the transition temperature for a = 0, i.e., 
kBTcQ = ecexp(-l/At) with 

Atd(k) =-7V(0)(Vuk'd(k'))k,. (17) 
In the following we assume q || z, for simplicity. Then 
gk+a is independent of q, gk+a = ^J\^{-ky,k^,Q). In 
the next section, we calculate the transition temperature 
as a function of a in q = for the order parameters 
shown in Table I which are the basis functions for the 
Civ point group suggested by A. Sergienko and S. H. 
Curnoe^'' and P. A. Frigeri et al.^ 



d(k)=Vl'2(1/kF)(-ky, kx, 0) 

d(k)= \/3(1/kF)(kz,0,0) , 

v^(1/kF)tO, kz, 0) 
d(k)-\/^(1/ki=)(kx, -ky, 0), 

V572(1/kF)(ky, k:<, 0) 

---- d(k)= (1/kF)(kx, ky, kz) 

d(k)= x/372(1/kF)(k:<, ky, 0) 

d(k)= ^/772t1/kF)tkJ, k^t, 0) 




ll ■ 

0.5 



3. Transition temperature of the uniform state 

Fig. 1 display the transition temperature as a func- 
tion of a for some types of spin triplet superconduc- 
tivity with q — Q. The transition temperature of the 
state with d(k) || gk does not depend on a. As shown 
in Appendix, we obtain that the transition temperature 
for d(k) = V3^(fc,,0,0) and d(k) = y3^(0,fc„0) 
(the two-dimensional representation E) approaches to 
zero as proportional to , for large a . The tran- 

sition temperatures for d(k) — yj^^{kx,—ky,0) and 

3 1 



d(k) = J ^ — {ky,kx,0) also go to zero as proportional 



to 



for large a. The transition temperature for 
d(k) = -^{kx, ky, kz) is zero at a « O.SlnkBTco and 

the transition temperature for d(k) = \J\-^{kx,ky,Q) 
is zero at a « 0.267rfcBTco. Interestingly, the transi- 
tion temperature for d(k) — ky,0) depends 
on a in the similar way as the transition tempera- 
ture for d(k) = fcy, 0), but it approaches 
to zero as ( , "^p )~'^^ for large a. Curves of T, vs a 
for d(k) = -^{kx,ky,kz), d(k) = ^^{kx,ky,0) and 

d(k) = \f^T^{kxj ky, 0) are reentrant and they are rem- 
iniscent to the transition temperature for the singlet su- 
perconductivity as a function of external magnetic field, 
where only the Pauli pair breaking effects are taken into 
account and the orbital effects are neglected.^''' 

Here we summarize the known results for the Pauli 
pair breaking effect in the spin singlet superconductiv- 



Fig. 1. Transition temperature as a function of the strength of 
the spin-orbit coupling a for some types of spin triplet super- 
conductivity with q = and = tj ^-j^{—ky, kx,0). When 
d(k) II gk, Tc is independent of a. For the states with d(k) oc 
(kz,0,0),(0,kz,0),{kx,-ky,0) and {ky,k^,0), Tc is reduced by 
a as Tc oc l/a. For d(k) oc {kx,ky,kz) and (kx,ky,0) Tc be- 
comes zero at the critical value of a. For d(k) oc {k^,ky,0), Tc 
is reduced sharply near a ss O.SnkgTcO and Tc oc a"^^ at large 
a. 



singlet superconductivity Fermi surfaces is split into two 
concentric spheres by the external magnetic field. If we 
assume the second order phase transition into the super- 
conducting state with q = 0, the reentrance is obtained. 
However, this reentrance cannot be realized, since the 
first order transition to the superconductivity with q = 
is expected in the region where the reentrance might ex- 
ist. The first order transition into the q = supercon- 
ductivity from the normal state does not realized neither, 
since it is covered by the second order transition into the 
pairing with finite q (FFLO state). 

The similar situation is expected in the case of the 
pair breaking due to an anisotropic spin-orbit coupling. 
In case of spin-orbit coupling electron energies for the 
Hamiltonian (eq. (4)) are 



e± =^k±ayfc^ + A:2, 
and the eigenstates are 



1 ( ' 

ip± = —j= ky-ik^ 

v2 V V^S+fef 



(18) 



(19) 



ity 



17, 19 



In the case of Pauli pair breaking for the spin where the upper and the lower components are the up 
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(a) 



kz 



kx=0 





( 






1 


/ J 


k 



Fig. 2. Fermi surfaces in the presence of the spin-orbit coupling. 
Spherical Fermi surface (^^ = 0) is split into two surfaces, = 
and e_ = 0. 



and down spins. The Fermi surface is split as shown in 
Fig. 3. In the present case Fermi surfaces is split due to 
spin-orbit coupling as two spheres with different centers 
as in Fig. 2 and Fig. 3. The first order transition into the 
spin triplet state with q = might occur, but we do not 
consider that possibility in this paper. The second order 
transition into the state with q 7^ in the next section. 

4. Nonuniform state of the spin triplet pairing 
due to anisotropic spin-orbit coupling 

In Fig. 4, we plot transition temperatures of the states 
d(k) — \J\-^{kx,ky,Q)^ as a function of a for some 
values of g, which is defined by 



-q- 



(20) 



For q = 0, the transition temperature shows a reen- 
trance near a/ (nkBTco) ~ 0.3 (oc < a < ao, ~ 
0.2637rfcBTco, and « 0.3207r/cB'7co)- The reentrance 
never occurs, since the transition curve is calculated by 
assuming the second order transition. When the temper- 
ature is lowered at fixed a (ac < a < ao)i the order 
parameter is finite below the higher transition tempera- 
ture and the lower transition temperature is spurious, as 
in the case of the Pauli pair breaking in the spin singlet 
superconductivity. For the finite value of q the critical 
value Qfc at Tco — becomes larger. However, the criti- 
cal value of a never exceed ao ~ 0.327rfcBT'co for g = 0. 
Therefore, contrary to the Pauli pair breaking in the spin 
singlet pairing, the spin orbit coupling will not stabilize 

the nonuniform state of d(k) = \J\Y^{kx^ ky,0). 

We obtain the similar a dependence of Tc for the su- 
perconducting state of d(k) = -^{kx, ky, kz) as shown 
in Fig. 5. The reentrance seen in the region 0.263 < 
, "rp < 0.360 when g = 0. When a 0, the reentrant 



(b) 



kz=0 



— »— ^ 



kx 



Fig. 3. Cross section of the Fermi surface with the plane 
(a) and the plane fcz = (b). 



region becomes small, but the superconducting state is 
not realized in the region a > 0.360. Then the nonuni- 
form superconducting state is not stabilized for this order 
parameter. 

For the superconducting state with d(k) = 
fcy, 0), we find that the nonuniform state is 
stabilized in the region , " > 0.312 as shown in Fig. 6. 
In this type of the superconductivity is reentrant when 
a < 0.312 and (x a~^^ for Tc ^ Tco and a ^ nkBTco- 
As seen in Fig. 6, the nonuniform state with < g < 4.5 
is stabilized in the when a > 0.312. 



5. Conclusion 

We study the pair breaking effect due to the spin- 
orbit coupling for the spin triplet superconductivity. 
We have studied various type of spin triplet supercon- 
ductivity, which is characterized by the d vector, the 
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dCk)= v^CI/kF) Ckx,ky,0) 




a/7ikBTco 



Fig. 4. Transition temperature as a function of a for triplet 
superconductivity with = \J \^{~^yt^x,0) and d(k) = 



states with d(k) 



^lyiO)? f^ikxiky^kz), and 



/c^, 0) show the reentrant curves of Tc as a 
function of a. Although the non-uniform state [q ^ 0, 
FFLO state) is not stabihzcd in the former two states, we 
find that the nonuniform state is stabihzcd in the latter 
state. The region of a where the nonuniform state is sta- 
bilized is rather small. In this study we have assumed the 
spherical Fermi surface. If we consider the Fermi surface 
with different shape with lower dimensionality (quasi- 
two-dimension or quasi-one dimension) or better nesting 
condition of the Fermi surface, FFLO state is expected 
to appear in wider region of parameters, as is the case in 
the ordinary FFLO state caused by the Pauli pair break- 
ing. 

If the FFLO state is realized in the absence of magnetic 
field, the amplitude of the energy gap varies in space and 
the specific heat and other properties will depend on tem- 
perature as a power low. The experimentally observed 
properties of CePtaSi, i.e. spin triplet superconductivity 
with power-low temperature dependences of penetration 
depth and thermal conductivity, are consistent with the 
nonuniform spin triplet state due to the spin-orbit cou- 
pling. 

Acknowledgement 

We thank H. Shimahara and N. Hayashi for useful dis- 
cussions. H.T and H.K. thank M. Sigrist, D. Agterberg 
and Y. Yanase for fruitful comments. 



d(k) = (1/kF) (kx, ky, kz) 




O 0.1 0.2 0.3 0.4 

a/itkBTDd 




a/;ikBTc 



Fig. 5. Transition temperature as a function of a for triplet 



superconductivity with 
J— (fci, ky, kz). 



'l^{-ky,k^,Q) and d(k) 



Appendix: Transition temperature in the case 
of strong pair breaking due to spin- 
orbit coupling 

In this appendix we calculate the transition tempera- 
ture of the uniform triplet superconductivity in eq. (14) 
for q = and ^ 0. For q = 0, eq. (14) becomes 



Tr 



ln(-^) = -2'KkBU[\ d(k) 1^ -(gk • d(k))^ 



rie-l 



i\ujn\ 2i I w„ I +2a \ gk 



1 



2i\ujn \ -2a I gk 



(AT) 



We set pk = -St^- As we assume that Fermi surface is 

Trks-lcO 

spherical, 

kx ~ kp sin 6 cos (p (A--2) 

ky — kp sin0 sin (/) (-A-3) 

kz = kF cos 6 (A-4) 

and the average over the Fermi surface is performed as 



dOsinO ■ 



(A-5) 



Then the transition temperature is obtained by 



T, 



x-Re 



hi '^^j, ^(^) I' -(gk • d(k))^ 



^K + ^)-*(i) 



(A-6) 
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0.3 

a/TiksTco 



Fig. 6. Transition temperature as a function of a for triplet 
superconductivity with gk = \J ^'i~(~ky, ^x,0) and d(k) = 



Substituting eq. (A-9) into eq. (A-8), we obtain 

T 1 f^^ 3 

ln(-^) = -— / dcj, dz-z\l + cos2<j>) 
Jco 47r Jo J-i ^ 



X In ( 



v/6(l-z2)g7Q,^ 



dz 



2, fV^e^a 1 2, /I 



We perform the integration over using 



dzz''ln(l - z^) = -(ln8-4), 
9 



and we obtain 
T 

H^) = 



ln( 



irknT, 



ln(2\/6, 



(A-IO) 
(A-11) 
(A-12) 



Finally we obtain that the transition temperature for 
Tc ^ TcO is obtained as 



(A-13) 



0.435 



cO 



(A.14) 



i.e. Tc is reduced as a ^ in the strong spin-orbit coupling 
case (a » 1). 



where -^{uc + \) is the digamma function. Since > 1, A. 2 d(k) = \ ^^iK, -ky, 0) 
we set + \) ~ *(rtc + | + i^) ~ in eq. (A-6). At 
low temperatures (Tc ^ Tco), pk 3> 1. Then 



In this subsection we study the case of d(k) 



Re(vI/(i+^^))«Re(*(^^) 



3 1 



Re 



-H^) (A-7) 

and the transition temperature in the case of Tc ^ Tco 
is obtained by 



-^{kx, ~ky,0). In this case we obtain 



I d(k) |2 -(gk-d(k))2 = -(l-cos2 0)(l+cos4</)), (A-IS) 

and we perform the integration over and k^ in eq. (A-8). 
We use 



and 



[ dz\n{l ~ z^) = -2 + ln4 
Jo 



dz7?h\(\~ 7?) = -(ln8-4). 



(A.16) 



(A-17) 



""^Tco^^ 4^ 



ln(2f 



27r 

I 

J -1 

a I gk 



d{cose) [I d(k) |2 -(gk • d(k))2]In this case we obtain, 

-1 

), (A-8) 



TrfcsTt 

where we have used ^"(1) = — ln(4e''') and 7 is the Euler 
constant. 

A.l d(k) = V3ji^(A:„0,0) 

First, we study the order parameter d(k) = 
^3-^(^2,0,0). In this case 

I d(k) |2 _(gk • d(k))2 = ^ cos2 e{l + COS20). (A-9) 



T 1 3 /■^'^ 



x(l + cos40)ln(^^iip5^) 



1, ,V6e 
-ln( 



2 ' TTfcijT, 

^.^(31n2-W). 



nksTc 



-)--(-2 + ln4) 



This equation is rewritten as 



Tco' VfcsTcV "^''^ 2Vf ^' 



(A-18) 



(A-19) 
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Finally we obtain 



1 TTfesTcO 



« 0.263^^^. 



(A-20) 



A.3 d(k) = y|-^(/c^,fc„0) 
For d(k) = y^j^(fcj/5fcx)0) we obtain 

I d(k) |2 -(gk-d(k))2 = ^(l-cos2 0)(l-cos#). (A-21) 

Since the term proportional to cos 40 vanish in the inte- 
gration over (j), we get the same a dependence of as 

the case of d(k) = ^J\-^{kx,—ky, 0). 

A. 4- d(k) = -^{kx, ky, kz) 
When d(k) = {kx, ky, k^), gk • d(k) = and 

|d(k)|2-(gk-d(k))2 = l. (A.22) 

We substitute eq. (A- 22) into eq. (A- 8), and obtain 



rri -j p2,7T p 1 



V'6(l - z^)e^c 



■) 



= -H^^fr) - ln(2N/6e^-i). (A.23) 



We obtain 



a 



0.31 



(A-24) 



from eq. (A-23), i.e. Tc = at a w O.SlTrfcsTco. 

^.5 d(k) = ^^(fc,,fc„0) 

When d(k) = ^J-{kx,ky,0), gk • d(k) = 0, 

I d(k) |2 -(gk • d(k))2 = ^(1 - cos^^). (A.25) 

We substitute this into eq. (A-8), 

, 3, V6(l-^2)e7a 



ln(;^) = --^ / / dzhnC- 
^T^o' 4nJo ^y_i 2 ^ 



= -H^^) - ln(2x/6e^-t). (A.26) 



Then we obta: 



1 



i.e. Tc = at a w 0.267rA;BTco. 

^.^ d(k) = Jf^(fc3,fc3^0) 



0.26, 



(A-27) 



When d(k) = ^lJ^{klk^y,0), 
|d(k) p-(gk.d(k))2 



= — (l-cos^6i)'^(19-|-12cos#-|-cos80). (A-28) 



Substitute it into eq. (A-8), we obtain 



""^T^' 47r64 



<.27r 



2n3 



dzil-k^) 



/6fl — z^'le'^'a 
X (19 + 12cos4(^ + cos80)ln( ^ , )■ 



irksTc 



(A-29) 



Using 

/■ 

Jo 



dkzil - z^f\n{l - z^) 



2552 16 



3675 35 



ln4, (A-30) 



we obtain from eq. (A-30) 

, , Tc , 19,, ,2y/6e^a,. 19 319 

^"(7^) = -20N^^)) + I^-^^ (A-31) 

which is written as 

Finally we obtain 

~ ^ 2x/6 ^ ^TrfcsTeo^ 



254!lMk') 



19 



5.08 X 10" 



12 / TrfceT'f 



cO 



19 



(A-33) 



The transition temperature is reduced as a when a ^ 
1. 
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